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ABSTRACT 



In a recent paper, Romero et al. (1999) have pre- 
sented the results of a positional correlation analysis 
of the low-latitude unidentified gamma-ray sources 
in the Third EGRET catalog with complete lists of 
galactic objects like supernova remnants, early type 
stars with strong stellar winds, and OB associations 
(considered as pulsar tracers). We have now carried 
out a study of those 42 sources at \b\ < 10° for which 
no counterpart was found at all. A variability anal- 
ysis shows that this sample contains a population 
of sources with high levels of variability at gamma- 
rays. The surface density of these variable sources is 
5 times higher than what is expected for unidentified 
AGNs seen through the Galactic plane. We discuss 
the origin of this presumed Galactic population of 
gamma-ray objects and the role that INTEGRAL 
could play in their physical identification. 
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the positional uncertainties, with Wolf-Rayet and Of 
stars (which are hot and massive stars endowed with 
strong supersonic stellar winds), 26 with OB associ- 
ations (usually considered as pulsar tracers) and 22 
with supernova remnants (SNRs). The probability of 
pure chance superposition with OB associations and 
SNRs is quite negligible (< 10 -5 ). For the stars, the 
case is suggestive although not so conclusive (proba- 
bilities < 10~ 2 ). Romero et al. (1999) also found 
a set of 42 gamma-ray sources for which there is 
not positional superposition with any known Galac- 
tic object capable to produce high-energy gamma- 
rays. The nature of this sub-group of sources remains 
a mystery. 

In this paper we present a variability analysis of the 
sample of unidentified 3EG sources near the Galac- 
tic plane that do not present positional correlation 
with potential Galactic counterparts. Such an anal- 
ysis can shed some light on the nature of the par- 
ent population, because radio-quiet pulsars and yet 
undiscovered SNRs are not expected to present high 
levels of variability on timescales of months. Based 
upon the results, we shall advance a working hypoth- 
esis on the parent objects. 



1. INTRODUCTION 



The Third EGRET (3EG) catalog of high-energy 
gamma-ray sources lists 271 point-like sources de- 
tected at E > 100 MeV (Hartman et al. 1999). 
Among these sources there are 66 high-confidence 
identifications with blazars, 7 confirmed gamma-ray 
pulsars, and about 170 unidentified sources. The dis- 
tribution of the unidentified sources with Galactic 
latitude shows a strong density gradient towards the 
Galactic plane, which indicates a strong contribu- 
tion from sources belonging to our Galaxy (see Fig- 
ure 1). Additionally, a group of mid-latitude sources 
has been recently found to be correlated with the 
Gould belt, a nearby region of active star formation 
(Grenier & Perrot 1999, Gehrels et al. 2000). 

The positional correlation between the 3EG sources 
in the Galactic plane and different types of Galac- 
tic objects was studied by Romero et al. (1999). 
They found that 10 sources are coincident, within 



2. VARIABILITY ANALYSIS 



There are 81 unidentified sources at |6| < 10° listed 
in the 3EG catalog. We shall focus our analysis on 
the sample of 40 sources with no known counterpart 

Zhang et al. (2000) have already presented a vari- 
ability analysis of the sources that are correlated 
with OB associations and SNRs, finding that most 
of these sources could be non- variable. These au- 
thors have shown that this group of sources could be 
mostly originated in a population of young and/or 

1 We have excluded 2 sources that are almost certainly ar- 
tifacts associated with the proximity of the very bright Vela 
pulsar (these sources do not show up in a map which excludes 
the Vela pulsation intervals) and hence the difference with the 
number of 42 sources mentioned by Romero et al. (1999). 
There are other 4 possible artifacts that are positionally cor- 
related with SNRs. 
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Figure 1. Distribution of unidentified gamma-ray 
sources with Galactic latitude, for the interval \b\ < 
10°. 



radio-quiet gamma-ray pulsars. We aim here at es- 
tablishing the variability properties of those sources 
not included in Zhang et al. (2000) analysis. 

In order to carry out the variability analysis we define 
a mean weighted value for the EGRET flux as: 
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Here, N vp is the number of viewing periods for each 
gamma-ray source. We take into account only single 
viewing periods. F(i) is the observed flux in the i- 
period, whereas e(i) is the corresponding error. For 
those observations in which the significance (\/TS in 
the EGRET catalog) is greater than 3a, we took the 
error as e(i) = F(i)/y/TS. However, many of the 
observations are in fact upper bounds on the flux, 
with significance below 2a. For these ones, we as- 
sume both F(i) and e(i) as half the value of the up- 
per bound. This is a rather conservative assumption 
that is not expected to result in an overestimate of 
the variability levels, as it could be the case if zero 
flux is directly adopted. We then define the fluctua- 
tion index [i as: 



/i = 100 x a sd x (FY 



(2) 



In this expression, <7 s d is the standard deviation of 
the flux measurements, taking into account the pre- 
vious considerations. 

In order to remove as far as possible spurious vari- 
ability introduced by the observing system, we com- 
puted the fluctuation index fi for the 7 confirmed 
gamma-ray pulsars (i.e. those listed in the 3EG cat- 
alog plus the sources 3EG 0634+0521 and 3EG 1048- 
5840, which were recently identified by Cusumano et 
al. 2000 and Kaspi et al. 2000, respectively). We 



adopt the physical criterion that pulsars are -i.e. by 
definition- non-variable gamma-ray sources. Then, 
any non-null /x-value for pulsars is attributed to ex- 
perimental uncertainty. We then define a statistical 
index of variability, 7, as 



<C jl > pulsars 



(3) 



Then, non-variable sources are defined as those for 
which I < 1 + la. Sources with 7 > 1 at a 3a 
level are classified as variable (la = 0.5). Sources 
with 7 > 1 at less than 3a are dubious cases and we 
cannot conclude about their variability within the 
present observational accuracy. 

We adopt this criterion, previously used in blazar 
variability analysis (e.g. Romero et al. 1994), instead 
of the similar one used by McLaughlin et al. (1996) 
in order to have a direct comparison with the spu- 
rious statistical variability shown by pulsars. Since 
the 7-index establishes how variable is a source with 
respect to the pulsar population, it can be considered 
as indicative of how confident we can be about the 
possible physical variation in the gamma emission of 
the sources. 



3. RESULTS 



The result of our analysis is presented in Figure 2 in 
the form of an histogram. We have fitted the distri- 
bution with a Gaussian curve in order to estimate the 
position of the maximum. The peak of the histogram 
occurs at 7 — 2.0, which means that there is an im- 
portant contribution from sources that display signif- 
icant variability over timescales of months. About 30 
% of the sample (12 sources) presents indices above 
2.5, which are clearly variable at more than a 3a con- 
fidence level. There are sources with indices as high 
as 7 = 8.9, like 3EG J1735-1500. There are also 18 
dubious sources (45 %) with 1.5 < 7 < 2.5. Only 
10 out of 40 sources (25 %) are clearly non-variable. 
The variability histogram for the sources in the sam- 
ple is very similar to that obtained for AGNs in the 
3EG catalog (see Figure 3), which is a well known 
variable population. 

Our results are, in general terms, similar to those 
recently presented by Tompkins (1999), who re- 
analyzed the EGRET data to take into account not 
only all sources included in the 3EG catalog, but 
also the 145 marginal detections not included in the 
final official list. Two of the sources pointed out by 
Tompkins as the most variable ones at low latitudes 
are included in our sample with indices of 3 and 5.3. 
Although the dubious cases trend to be non variable 
in the most comprehensive Tompkins analysis, our 
method is, we think, still the best approach when 
data from the published catalog are used. 

In Figure 4 we present a plot of the variability ver- 
sus the photon spectral index T (defined such that 
F(E) oc E~ r ) for sources in our sample. In order 
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Figure 2. Distribution with variability index for 
sources that are not correlated with known Galactic 
objects. 
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Figure 3. Distribution with variability index for 
sources that are identified with AGNs in the 3EG 
catalog. 



to provide a comparison we have also included in 
this plot the population of known gamma-ray pul- 
sars and the population of identified AGNs. All 
known gamma-ray pulsars are located between the 
solid horizontal lines. It is interesting to note that 
those unidentified sources with the steepest spectra 
seem to have the highest variability indices. The 
probability of this effect being the effect of chance is 
~ 8 %, not too low, but sufficient as to encourage 
further studies. 

Figure 5 shows the histogram with distribution of 
the high-energy photon index T for the sources of our 
sample. Notice that known gamma-ray pulsars have 
hard spectra at E > 100 MeV, with indices typically 
r < 2.15 (e.g. Thompson et al. 1994), contrary 
to most of the sources here under discussion, which 
present an average index <T>=2.4±0.2. 
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Figure 4- Variability vs. photon spectral index for 
the sources under study. Known gamma-ray pulsars 
and identified AGNs are also shown for comparison. 



The nature of the variable population of unidentified 
sources is a very intriguing problem that might be 
addressed by the INTEGRAL mission. 



4. MAGNETIZED BLACK HOLES 



The variability analysis presented in the previous sec- 
tion clearly shows the existence of a galactic popu- 
lations of variable gamma-ray sources that do not 
show up significant emission at other wavelenghts. 
The surface density of variable sources in the sam- 
ple here under study is five times higher than what 
is expected for unidentified AGNs seen through the 
Galactic plane. 

Recently, Punsly et al. (2000) proposed that variable 
galactic sources, which from its gamma-ray proper- 
ties resemble very much AGNs, could be isolated 



Kerr-Newman black holes. Punsly (1998) has shown 
that a magnetospheric charge may be supported in a 
disk or ring around a charged, rotating black hole of 
a few solar masses. The total configuration is stable 
and has zero net charge, hence the hole should not 
discharge quickly through accretion from the diffuse 
ISM. Unlike the case of pulsars, magnetized black 
holes do not present solid surfaces so they do not 
produce thermal X-ray emission. Strong magnetized 
bipolar winds are expected in the form of collimated 
electron-positron jets (Punsly 1998), where gamma- 
ray emission is generated through inverse Compton 
scattering of synchrotron photons and pair annihila- 
tions. 

Detailed calculations of the spectral energy distri- 
bution of these objects (Punsly et al. 2000, Combi 
et al. in these proceedings) show that at MeV en- 
ergies the annihilation luminosity exceeds the self- 
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Figure 5. Histogram with the distribution of the high- 
energy photon index for the sources included in the 
present study. Notice that most sources have steeper 
indices than those presented by known gamma-ray 
pulsars. 



~ 10 35 erg s -1 (Combi et al., these proceedings), so 
if it occurs at, say, 4 MeV in a source at ~ 1 kpc, 
it should be clearly detected (see Combi et al., these 
proceedings). The exact position of the peak in the 
spectral energy distribution can be used to infer the 
Doppler factor of the jets, whereas its measured in- 
tensity can be applied to the calculation of the pair 
density in the relativistic flow. The annihilation lu- 
minosity is given by (Punsly ct al. 2000): 



t ann 
L 



= (3/32)cr T c(r - iyn A m e c A V x 



+ 



2(2r- 1) 



(r- 1/2) (r + 1/2) r 2 (r-i) 



(4) 



where ut is the Thomson cross section, V is the vol- 
ume where the annihilations occur, n is as before 
the number density of electron-positron pairs, T is 
the high-energy photon index, and m e is the elec- 
tron rest mass. 

Frequent observations of these sources with IBIS 
will be useful to establish the short-term variability, 
which in turn can shed light on the involved plasma 
instabilities. High-resolution radio observations to- 
wards those sources detected by IBIS could find the 
synchrotron radio signature of the jets, that should 
appear as a point-like non-thermal source of a few 
tens of mJy at 5 GHz (Punsly et al. 2000). 



Compton luminosity, producing a break in the spec- 
trum and a steepening at high-energies. The sources 
with the steepest gamma-ray spectra are those whose 
jets have the highest pair number densities. Dense 
jets, in turn, are prone to undergo plasma instabil- 
ities. Synchrotron cooling will most likely cause ra- 
diative losses of the random transverse (thermal) en- 
ergy in a magnetized jet. One of the instabilities 
resulting from the expected pressure anisotropy in a 
jet is associated with the shear Alfven wave: the fire- 
hose instability. The fastest growing modes of this 
instability have a growth rate cx n 1 / 2 , where n is the 
pair density. Hence, very dense jets with large an- 
nihilation luminosities and a steep spectral indices 
tend to be very firehose unstable. The wobble im- 
plies variations of the Doppler factor of the beamed 
emission, 8, and consequently of the luminosity of the 
jet, which has a dependency L cx 5 3 . Hence, those 
sources with the steepest gamma-ray indices would 
display the highest variability indices, as observed in 
our sample of unidentified gamma-ray sources with 
no positional counterpart at lower frequencies. 



5. ADDITIONAL COMMENTS 



The magnetized black hole hypothesis for the par- 
ent population of the low-latitude, variable gamma- 
ray sources can be tested by the IBIS instrument of 
INTEGRAL mission. The spectra of these objects 
should exhibit a broad peak at a few MeV, similar 
to that observed in MeV blazars (e.g. Bloemen ct 
al. 1995). The expected luminosity of this peak is 



6. CONCLUSIONS 



In this paper we have estimated the variability of 
a sample of 40 low-latitude, unidentified gamma- 
ray sources that are not positionally correlated 
with known potential Galactic emitters. We have 
shown that within this set there exists a group of 
sources which display high levels of variability over 
timescales of months to years at gamma-rays. These 
sources are unlikely to be pulsars or undiscovered 
SNRs. The INTEGRAL satellite will provide a 
powerful tool through its IBIS imager to probe the 
nature of this population of energetic objects and 
test whether the hypothesis here presented, namely 
that the gamma-ray emission is produced by isolated 
magnetized black holes, is correct or not. 
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